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At the Fermilab Tevatron, final states with three tagged 6-jets 
could play an important role in searches for a Higgs boson 
with mass in the range 100-300 GeV. These signals arise from 
gb fusion and we demonstrate their observability in the limit 
of a large 6-quark Yukawa coupling. Rather promising dis- 
covery limits on such a coupling are obtained and consequent 
effects on the parameter space of the Higgs-boson sector in 
the MSSM are discussed. 

PACS number(s): 12.60.-i, 14.80. Cp, 14.65.Fy, 12.38.Bx 

The Higgs boson is the last ingredient of the Stan- 
dard Model (SM) that still awaits discovery. It is also a 
crucial component since the mechanism for electroweak 
symmetry-breaking — and hence the masses of weak 
gauge bosons and fcrmions — arises from its interac- 
tions. There exist strong indications from the fitting of 
electroweak parameters to available data that the Higgs 
boson has a relatively small mass jj]], which makes it 
likely to be seen either at the Fermilab Tevatron (in its 
future runs) or at the CERN LHC (when it becomes op- 
erational) . Curiously, although failure to find a Higgs bo- 
son at all is likely to require serious rethinking of current 
theories of electroweak interactions, establishing its exis- 
tence could also raise a rash of questions which might not 
be easy to answer. For example, one could ask whether 
a discovered Higgs boson is the one predicted in the SM 
or whether it is an ingredient of some extended theory 
with extra symmetries. The easiest way to settle this 
question would be by the detection of (an) additional 
Higgs boson(s) such as is (are) predicted in several of 
these models. Alternatively, on discovering a Higgs bo- 
son, we could determine its couplings and compare these 
with theoretical predictions within a given model. With 
the higher integrated luminosity expected in Run II of 
the Tevatron, the requisite precision in the measurements 
should be available, at least for favourable regions of the 
parameter space. A truly dramatic way of discovering 
new physics beyond the SM would be the observation of 
Higgs-boson signals where none are predicted in the SM. 

Enhanced Yukawa couplings for the 6-quark in some 
models beyond the SM open up many such interest- 
ing new channels. At the CERN LEP, for example, 
e + e~ — > Z* — > bbA could help [g] extend the Higgs-boson 
discovery limits beyond the expectation from the study of 
Higgs-strahlung processes, which are the main avenue of 
Higgs-boson searches at present. At the Fermilab Teva- 
tron, similar processes have been studied, some recent 
efforts being those of Refs. pl-pl]. The consensus appears 



to be that detection of the SM Higgs boson radiated off 
a heavy quark would be difficult, but (non-SM) Higgs 
bosons with enhanced Yukawa couplings would be much 
better candidates for a study in 46 and bbrr final states. 
In this letter, we propose the related process 
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as a possible production channel for a generic Higgs bo- 
son, 4>. While the analogous interaction for the charged 
Higgs boson (in two Higgs-doublet models) has been 
studied extensively in the context of the CERN LHC, 
this particular process has been neglected so far, proba- 
bly under the impression that the cross section would be 
hopelessly small. This is, in fact, true for the SM Higgs 
boson, unless we consider the TeV-33 option at the Fer- 
milab Tevatron. However, we shall show that it need 
not be the case - even in Run II - for theories with an 
enhanced (j)bb coupling. 

The Yukawa interaction can be written C^b = hb4>bb, 
or C^ b i — ihb^bj^b, depending on the CP assignments of 
the scalar <fi. Given such an interaction, the cross-sections 
(to leading order) for the processes in Eq. ([!]) scale as 
the ratio (hb/hf M ) 2 , irrespective of the CP properties 
of the field (p. The precise value of the ratio hb/hf M 
depends, of course, on the model. In the minimal su- 
persymmetric extension of the SM (MSSM), for exam- 
ple, hb/hf M — tan/3, — sincv/ cos/3 and cosa/cos/3 for 
the pseudoscalar (A ), light scalar (h°) and heavy scalar 
(H ) Higgs bosons respectively, where the angles a,j3 
have their usual meanings Q . In the large tan (3 limit, the 
coupling to the A (obviously) and the H° are strongly 
enhanced; that to the h undergoes a more modest in- 
crease. Similar enhancements of the 6-quark Yukawa cou- 
plings can occur in composite models such as those with 
topcolour-assisted technicolour . 

Once a <j) has been produced, various decay chan- 
nels are available to it. Again, there is considerable 
model dependence. For the SM H°, the major two- 
body decay modes are H° —> bb,cc,r + T~ with relative 
branching ratios of 23 : 2 : 4 approximately M. While 
these modes dominate for low Higgs-boson masses, the 
H° — » WW* mode becomes increasingly important for 
m H > 100 GeV. In the MSSM, on the other hand, the 
WW* (or ZZ*) mode is inaccessible to the A . Con- 
sequently, Br(A° — » bb) <; 90% (unless ^4° can decay 
into a pair of superpartners). Similarly, in the event of 
tan/3 3> 1, there are just three important decay modes 
for the scalar h°, namely h° — > bb, t + t~ ,WW* , the first 
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always accounting for more than 70%. While the WW* 
mode may also be profitably used as a probe of Higgs- 
boson interactions, we prefer to concentrate on the sim- 
pler two-body decays of the Higgs boson into fermions. 
Thus - for the SM and MSSM at least - the final states 
of greatest interest in gb fusion are 6t + t _ and 666 |9j. Of 
these two, the former has smaller backgrounds, whether 
from QCD or from associated Z-production. However, 
Br((f> — > t + t~) is suppressed approximately by a factor 
of 6 or more over most of the parameter space. More 
importantly, invariant mass reconstruction for a r-pair is 
difficult because of multiple neutrinos arising in tau-pair 
decays. As mass reconstruction turns out to be a major 
tool in the isolation of a Higgs-boson signal, we shall not 
comment on the t + t~ channel any further |]l0| . 

Thus, our choice of the final state is (bbb) or (666), 
which we generically denote by 36. The signal cross sec- 
tion may be written as a(gb — > 36) = lZ 2 asM(gb — > 36), 
where 
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The advantage of using 1Z as a free parameter is obvious - 
it contains the entire model-dependence of the cross sec- 
tion and hence enables us to make a model-independent 
study of the 36 signal. 

At this point, it becomes necessary to comment on the 
size of hf M . The low-energy value can be inferred from 
the pole mass, for which we use m& = 4.3 GeV [0. At 
large momentum transfers, QCD corrections can be im- 
portant. Since the complete set of corrections have not 
been calculated, we include only the suppression due to 
the running of TO&. As additional QCD corrections usu- 
ally tend to increase the cross-section, our results should 
be regarded as a conservative estimate of the signal. Our 
expressions are consistent with those in Ref. G2|. 

The backgrounds to (fy) arise from two main sources: 
(i) 'authentic' 36 events from QCD and/or weak inter- 
actions; (m) spurious events of the type 26 + J, 6 + 2 J 
or 3 J, where J denotes a non-6 jet misidentified as a 6 
jet. Herein lies the advantage of considering a 36 sig- 
nal: misidentification probabilities are usually low and 
'authentic' 36 backgrounds carry the same suppression 
from the 6-quark flux as the signal. In contrast to this, 
46 backgrounds H could be generated by QCD processes 
arising from valence quarks or gluons, which have enor- 
mous fluxes by comparison. 

The large number of diagrams contributing to the 
background are calculated using the helicity amplitude 
package MADGRAPH |U5). To estimate the number of 
events and their distribution(s), we use a parton-level 
Monte-Carlo event generator. For the parton densities 
in the proton, we use the CTEQ3-M structure functions 
as incorporated in the package pdflib |L4j . 

Since the QCD background is populated mostly at low 
transverse momenta (pr) and high rapidity (r/) of the 



jets, we demand that the final state be composed of ex- 
actly three hard jets (J) with p 3 T > 20 GeV, \r)j\ < 2. 
We also require that the angular separation of the jets 
be substantial, i.e. ARjj = ^/(A?^-.,) 2 + (A<fijj) 2 > 1.0 , 
adapting the well-known cone algorithm for jet separa- 
tion to a parton-level analysis. While ARjj > 0.7 is usu- 
ally considered sufficient for jets to be separable, the more 
stringent cut (especially for the two softer jets) eliminates 
a significant portion of the QCD background without af- 
fecting the signal at all. To eliminate the background 
from gb — > bZ° , we demand that all events where the in- 
variant mass rriij satisfies 80 GeV < to^ < 100 GeV, 
for any of the three possible pairings (ij), should be 
rejected from the analysis. Similarly, a requirement of 
TOjj > 10 GeV helps to further reduce the QCD back- 
ground. Once these kinematic cuts are applied, we are 
in a position to utilize a particular feature of the sig- 
nal event topology. The final-state 6 in the production 
process gb — > b(j) tends to be collinear with the initial 6 
quark and hence (usually) to have a transverse momen- 
tum smaller than those of the 6's arising from the scalar 
decay cf) — > 66. If we label the 6-jets according to their 
Pt, thus: priji) > PtU?) > Prija) — the invariant mass 
to 12 of the pair of hardest jets will reconstruct to the 
Higgs-boson mass for a majority of the signal events. Of 
course, there will always be some events where the pair 
(12) does not originate from the scalar, but such configu- 
rations are subdominant and become progressively so for 
larger m^. Thus the TO12 distribution for the signal will 
exhibit a characteristic peak, illustrated in Fig. |l|. With 
the elimination of the Z events, the background does not 
show any such peaking. 
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FIG. 1. Distribution in the invariant mass of the pair of 
6-jets with the largest pr- Solid, dashed and dotted lines 
correspond to three different masses (marked) of the Higgs 
boson. 



We can thus use TO12 as a discriminator for Higgs- 
boson resonances. When looking for a <f> with a given 
mass, it is necessary to select only those events that 
lead to an TO12 close to to^. Closeness is usually quan- 
tified by the invariant mass resolution which is taken 
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to be |U max{10 GeV, ^0.8 GeV x m } for the entire 
range = 100-300 GeV. We note in passing that sim- 
ilar - but much weaker - peaks can be seen in the other 
two mass distributions TO13, 77123, although we have not 
exhibited them. 

Even with this kinematic selection procedure, the 
background is orders of magnitude larger than the signal. 
We demand, therefore, that all three 6's be tagged. Of 
course, progressively smaller fractions of 26+ J, 6+2 J and 
3 J states, respectively, will be misidentified as 36 states. 
To quantify this, for each individual jet, we use p5] a 6- 
tagging efficiency = 0.6 and misidentffication probabil- 
ity -Pmis = 0.005. This eliminates most of the (26+ J) and 
almost all of the 6 + 2 J and 3 J backgrounds while sup- 
pressing the signal (and the 'authentic' 36 backgrounds) 
by a factor of (0.6) 3 ~ 0.2. 

Using all these cuts and selection criteria, we are now 
in a position to compare signal with background. Fig. || 
shows the distribution in invariant mass m\2, assum- 
ing uniform bins of 10 GeV (l(| in the entire range 
mi2 = 100-300 GeV, for an integrated luminosity of 
L = 1 flV 1 . To illustrate the usefulness of this channel, 
we have considered a somewhat optimistic value 1Z = 50, 
for which it is obvious that the signal is significantly 
larger than the background. However, since it is only 
necessary for the signal to exceed the fluctuation in the 
background at (say) 95% confidence level (CL), the ac- 
tual reach in 1Z encompasses much lower values. 
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FIG. 2. The number of events in each 77112 bin for signal 
(points) and the total background (solid line). The individual 
contributions to the background are also shown (dashed and 
dotted lines). 

Fig. H also shows the individual channels contribut- 
ing to the background. For relatively low Higgs-boson 
masses, the largest background arises from 'authentic' 36 
processes, single misidentified non-6 jets come in second 
and double misidentification being strongly suppressed. 
Triple misidentification is severely suppressed and need 
not concern us further. For higher masses of the Higgs 
boson, however, the single-misidentified-jet background 
catches up with the 'authentic' one, so that the impor- 



tance of this channel hardly needs to be emphasized. 

We are now in a position to speculate on the mass 
reach of this mode. Concentrating on the appropriate 
mi2 bin, we can evaluate the minimum value of 1Z that 
would allow us to establish/exclude a Higgs boson with 
the corresponding mass. This is done in a straightforward 
manner by comparing the signal size (in that bin) with 
the statistical fluctuation in the background (in that bin) , 
calculated using Poisson statistics. 

In Fig. ||, we exhibit 95% CL discovery limits for three 
expected luminosities at the Fermilab Tevatron. If no 
excess in TO12 is seen, the parameter space above the 
curves can be ruled out at the appropriate confidence 
level. That the minimal value of 1Z grows with is ex- 
pected, since larger ra^ leads to small cross sections for 
the process (|l|). The minimum at ~ 110 GeV is an 
artefact of the kinematic cut designed to eliminate the 
Z background (in the absence of such a cut, of course, 
all the discovery limits would be much worse) . It is pos- 
sible, however, that a more sophisticated analysis could 
improve the bounds for ~ mz- The curve for a lu- 
minosity of 100 pb _1 represents the region which can be 
ruled out by existing data. 
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FIG. 3. Model-independent discovery limits for a given lu- 
minosity. The area above the curves can be excluded at 95% 
CL. 

In obtaining Fig. we have assumed that only one 4> 
lies in the mass bin under consideration. If this is not the 
case, or if there is more than one Higgs-boson state in the 
region under consideration, then the constraints originat- 
ing from all such </>'s will have to be compounded appro- 
priately. This could be parametrized as an enhancement 
in the effective value of 7Z, with the actual increase de- 
termined by the model parameters. The curves in Fig. ^ 
are thus conservative and, as we have already stressed, 
model-independent. 

In view of the intrinsic interest in the MSSM, we now 
pass from the general to the particular and derive ex- 
plicit constraints on this scenario. Before presenting our 
results, some preliminary remarks are in order. Given 
the current bounds on sfermions, it is clear that it is not 
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possible for a 100-300 GeV Higgs boson to decay into a 
pair of these. However, bounds on charginos and neu- 
tralinos are much weaker, and hence, we need to make 
some assumptions about the Wino mass parameter M2 
and the higgsino mixing parameter fi. In our analysis we 
choose two illustrative sets of such parameters. 

With the respective values of 1Z determined by rriA, 
tan/3, A/2, fJ-, the sfcrmion masses and the stop mixing 
parameter A t |12|, the three individual constraints can 
now be combined. The resultant bounds in the m^-tan (3 
plane are presented in Fig. Solid curves correspond 
to the case where all the superpartners are very heavy, 
thereby reducing the scenario to a constrained two-Higgs 
doublet model; dashed curves correspond to relatively 
light superpartners. With three non-trivial contributions 
(more than) offsetting any suppression due to branch- 
ing ratios, the bounds on tan/3 are analogous to those 
of Fig. U Interestingly, the dependence on the MSSM 
parameters is rather weak, as evidenced by the small dif- 
ference in the two sets. This is not unexpected as light 
Higgs bosons hardly decay into the superpartners, even 
when the mass parameters are set as low as 150 GeV. Ra- 
diative corrections in the Higgs sector, on the other hand, 
can significantly alter the individual couplings. On sum- 
ming over the three contributions, though, the residual 
effects are small and are not easily visible on the scale of 
Fig. ||. Our predictions are thus quite robust. 
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FIG. 4. 95% C.L. exclusion contours within the 
MSSM. Solid (dashed) lines correspond to the case 
mj = \fi\ = M 2 = 1 TeV (150 GeV). 

We may mention in passing that there do exist con- 
straints on the myt-tan/3 plane from direct searches at 
both the LEP collider and the Fermilab Tevatron |^-||] 
other than the 100 pb _1 curve shown here. However, 
most of these are outside the scale of the present fig- 
ure. It is obvious, therefore, that a 36 signal could be 
remarkably effective, if not in finding a Higgs boson of 
the MSSM, at least in putting stringent constraints on 
the parameter space. 

To summarize, then, we have discussed Higgs reso- 
nances in final states with three tagged 6-jets at the Fer- 



milab Tevatron. For Yukawa couplings of a Higgs boson 
to bb pairs - which are enhanced with respect to the 
SM coupling, we predict large signals in bins of invariant 
mass constructed using the hardest pair of jets. These 
are used to define a model-independent constraint on the 
parameter space, which is then translated to the case of 
a specific model - the MSSM. We show that this signal 
could lead to constraints on the MSSM parameter space 
that better present constraints by a significant margin. 
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